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We calculate the chirally improved (CI) Landau gauge quark propagator in two flavor lattice 
Quantum Chromodynamics and study thereon the effects of Dirac low-mode removal. The appli- 
cation of tree-level improvement of the propagator and tree-level correction of the lattice dressing 
functions removes the leading discretization artifacts. We find the dynamically generated mass in 
the infrared domain of the mass function to dissolve continuously with the reduction level and fur- 
thermore we observe strong suppression of the wave-function renormalization function for small 
momenta. 
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1. Motivation 

The lowest eigenmodes of the Dirac operator are related to the dynamical breaking of the 
chiral symmetry in QCD via the Banks-Casher relation. In [1, 2] we investigated the effects of the 
removal of the lowest Dirac eigenmodes from the valence quark sector on the hadron mass spectrum 
in a dynamical lattice QCD simulation. It was found that all hadrons, except for the pions, survive 
the aforementioned truncation. The masses of hadrons that can be transformed into each other by a 
chiral rotation became degenerate and surprisingly the value of these masses remained rather large: 
see for example Fig. 1 where we show the evolution of the vector and axial vector meson masses 
under Dirac low-mode removal. 

The removal of the lowest Dirac eigenmodes, i.e., the removal of the chiral condensate is 
expected to result in the loss of dynamically generated mass in the valence quarks. Thus, the naive 
expectation suggests rather light hadrons once the chiral symmetry has been restored. The study 
[3] aimed at shedding some light on the issue of mass generation in the context of Dirac low- 
mode removal: here we explored the valence quark propagators that have been used in the previous 
studies, in a gauge fixed setting in momentum space. Then it is possible to extract at each Dirac 
low-mode truncation level the mass function of the propagator which manifestly shows the amount 
of dynamically generated mass in the quarks. 
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Figure 1: The masses of the vector meson p and axial vector meson a\ under Dirac low-mode removal: the 
would be chiral partners become degenerate once the chiral condensate has been removed. 



2. Dirac low-mode truncation 

The lowest eigenmode of the Dirac operator are known to play a crucial role in the dynamical 
breaking of the chiral symmetry as stated by the Banks-Casher relation. Consequently, when re- 
moving the Dirac eigenmodes near the origin, the chiral condensate disappears and chiral symmetry 
becomes "unbroken" [1]. 

Consider the Hermitian Dirac operator D5 = 75D. We denote its (real) eigenvalues with fa and 
the corresponding eigenvectors with |v,). We write the Dirac operator D in terms of the spectral 
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representation of D=, : 

N 

D = £m^|v/)(v/|. (2.1) 
Then we can split the quark propagator S = D 1 into a low-mode part (lm) and a reduced part (red), 
S = £ Hf 1 \vi) (vi I 75 + £ \ii 1 |v,-) (vi | 75 = 5 lmW + S ied{k) (2.2) 

i<k i>k 

and hence we obtain the reduced part of the propagator by simply subtracting the low-mode part 
from the full propagator 

^red(A) = S — 5 lm (^ . (2.3) 

It is this reduced part of the propagator which we calculate within a gauge fixed setting in order to 
extract its momentum space dressing functions. 

We use dynamical two flavor gauge field configurations [4, 5] for our study where for the 
dynamical quarks as well as for the valence quark propagator that we investigate, the so called 
chirally improved (CI) Dirac operator [6, 7] has been adopted. We gauge fixed the configurations 
to Landau gauge with the cuLGT code [8] for lattice gauge fixing on graphic processing units 
(GPUs) with CUDA. 



3. The CI quark propagator from the lattice 

The continuum quark propagator at tree-level reads 

S(°)(p) = (/> + m)- 1 (3.1) 

where m denotes the bare quark mass. In manifestly covariant gauges the interacting renormalized 
quark propagator S(pi; p) can be decomposed into a Dirac scalar and a Dirac vector part 

S(n;p) = {irfA(n;p 2 )+B(ii; P 2 )y X (3.2) 

or equivalently 

S{pL;p)=Z{pL;p 2 ){i]/>+M{p 1 )y l . (3.3) 

Here we introduced the wave-function renormalization function Z{jl\p ) = \/A(p:;p ) and the 
mass function M(p 2 ) = B(pi;p 2 )/A(p:;p 2 ). 

The regularized lattice quark propagator St(p;a) can be renormalized at the renormalization 
point /I with the quark wave-function renormalization constant Z2(/i;a), 

S L (p;a)=Z 2 (lx;a)S(n;p). (3.4) 

The mass function M(p 2 ) is independent of the renormalization point pi whereas the wave- 
function renormalization function Z(fi;p 2 ) differs at different scales although can be related from 
different scales by multiplication with a constant. 

We extract the functions M(p 2 ) and Z(p 2 ) of the CI quark propagator from a the lattice. For 
details of the extraction see Ref. [9]. 
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The lattice quark propagator at tree-level 

sf\p)=(i4 + aMf\p)y l (3.5) 

differs from the continuum case due to discretization artifacts. The dressing function (p) equals 
one, by construction, and hence (p) equals at tree-level the mass function (p). 

We extract the CI lattice momenta ak{p) from the tree-level propagator and compare it to 
the analytically derived expression for the Dq\ momenta given in [3]. The lattice calculation is 
consistent with the analytical result, see the r. h. s. plot of Fig. 2. 
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Figure 2: The tree-level mass function (left in red) and lattice CI momenta (right) together with the corre- 
sponding analytical curves and the uncorrected interacting mass function (left in blue). 

The tree-level mass function aM^\p) which in the continuum is simply the bare mass m, is 
shown in the 1. h. s. plot of the same Fig. 2 (in red), again consistency with the corresponding 
analytical expression is manifest. As can be seen from the plot, aM^ (p) has a zero-crossing with 
aM[ 0) (0) ~ -0.333. 

When switching from the tree-level case to the non-perturbative full interacting case, the mass 
function aM^\p) of the 1. h. s. plot of Fig. 2 (red) turns into the interacting lattice mass function 
aM L (p) (blue). 

The shape of aMi(p) is similar to aMf\p) and is clearly altered by discretization errors. To 
get a handle on the leading lattice artifacts we perform improvement and tree-level correction. For 
the application of the Symanzik improvement program [10] to reduce the errors of the fermionic 
action to &(a 2 ) we refer to [3] and references therein. 

In order to blank out the lattice artifacts which are already present at tree-level, we aim at 
extracting the deviation of the interacting propagator from its tree-level form. To this end we 
simply divide the renormalization function Zi{p) by its tree-level form 

UP)^^\. (3.6) 

Z L (P) 

In order to apply an equivalent tree-level correction to the mass function of the form 

aM L ( P ) -+ (3.7) 
Mf\p) 
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we have to carry out an additive mass renormalization of the tree-level function B^\p) in order to 
avoid divergences because of the zero crossing of the tree-level function, 



aBi\p) ->■ aB^\p) + am a da 
where am^d is chosen such that B^' (0) = m, like in the continuum, i.e. 

am^dd = am — aB^ (0) « 0.344. 
Thus, the multiplicative tree-level correction for the mass function is 



aM L {p) 



amM L {p)A^(p) 

B L (P) + m add 



(3.8) 



(3.9) 



(3.10) 



4. Results 

We extract the quark wave-function renormalization function Zi{p) and the quark mass func- 
tion Mi(p) from the reduced propagators of (2.3) while varying the Dirac low-mode reduction 
level k form 2 to 512. We tree-level improve the modified propagators and apply the aforemen- 
tioned tree-level correction scheme. The dressing functions from the full propagators together with 
the dressing functions from the most extreme reduction (512 eigenmodes subtracted) are presented 
in Fig. 3. 
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Figure 3: Left: The corrected mass function from the original theory and after having subtracted the lowest 
512 Dirac eigenmodes. Right: The same comparison for the wave-function renormalization function of the 
quark propagator. 

The plot reveals amplification of infrared suppression of Zl(p) when subtracting Dirac low- 
modes while the range from medium to high momenta becomes not affected at all. Consequently, 
far traveling quarks are suppressed in this framework. 

The mass function Mi{p) gets strongly suppressed in the infrared when removing the lowest 
eigenmodes, the dynamic generation of mass ceases. In Fig. 4 we plot the deep infrared mass of 
the CI quark propagator from M^p^), at the smallest available momentum p m i n = 0.13GeV, as 
a function of the Dirac low-mode reduction level. We see that the value of the bare quark mass 
(horizontal line) is approximately reached at truncation level k = 256. Note that the reduction level 
k can be translated to an energy scale which is shown in the lower axis. 
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5. Conclusions 

We found that the dynamically generated mass in the valence quark propagator vanishes after 
a critical number of the lowest Dirac eigenmodes have been excluded from the simulation. This 
is in accordance with the Banks-Casher relation which relates the Dirac low-modes to the chiral 
condensate. Thus, the rather large mass of hadrons which are built out of such truncated quark 
propagators, see Fig. 1 and References [1, 2] cannot be explained by the dynamical breaking of 
the chiral symmetry. Interestingly, however, we find an intermediate Dirac low-mode truncation 
regime where the dynamically broken chiral symmetry seems to be restored, as indicated by the 
degeneracy of the masses of the chiral partner mesons p and a\ but, at the same time, the quark 
propagator still shows a significant amount of dynamically generated mass. Similar findings were 
obtained by the authors of [1 1] who performed a hadron spectroscopy on gauge field configurations 
from which center vortices have been cut out. 
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